Uniaxial and biaxial soft deformations of nematic elastomers 
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We give a geometric interpretation of the soft elastic deformation modes of nematic elastomers, with explicit 
examples, for both uniaxial and biaxial nematic order. We show the importance of body rotations in this non- 
classical elasticity and how the invariance under rotations of the reference and target states gives soft elasticity 
(the Golubovic and Lubensky theorem). The role of rotations makes the Polar Decomposition Theorem vital 
for decomposing general deformations into body rotations and symmetric strains. The role of the square roots 
of tensors is discussed in this context and that of finding explicit forms for soft deformations (the approach of 
Olmsted). 

PACS numbers: 61.41. +e Polymers, elastomers, and plastics, 61.20.Vx Polymer liquid crystals, 62.20.Dc Elasticity, elastic 
constants 



I. INTRODUCTION 

Nematic elastomers display three unique and related phe- 
nomena not found in conventional elasticity - large sponta- 
neous deformations, very large optical-mechanical response 
and soft elasticity. The last, shape change with little or no 
energy cost, is the subject of this paper. We give a geomet- 
ric interpretation of the soft modes described using fractional 
powers of tensors. We then discuss the soft modes of biaxial 
nematic elastomers. Since in contrast to conventional elastic 
solids, rotations play an essential role in the elasticity of ne- 
matic rubber, we conclude by discussing the related questions 
of breaking finite strains into symmetric shears and rotations, 
the Polar Decomposition Theorem for tensors, and the nature 
of square roots of tensors. 

Nematic elastomers have an internal, orientational degree 
of freedom in addition to those of ordinary elastic bodies. The 
anisotropy of molecular orientation induces shape anisotropy 
in the polymers that make up the elastomer. The switching on 
and off of this molecular shape anisotropy, either by temper- 
ature change[Jl|] or by illumination^, causes large (« 400%) 
mechanical shape changes. 

Rotation of this anisotropy, by imposed mechanical strains, 
is the extreme opto-mechanical effect that is observed. When 
accompanied by subsidiary shears and contractions which ac- 
commodate the changes of molecular distributions, the cost of 
the originally imposed strain is rendered to zero . This is what 
we call soft elasticity [^J. 



II. SOFTNESS IN LINEAR CONTINUA 

One can explore the related ideas of anisotropy, rotation 
and soft elasticity initially for small deformations and rota- 
tions, that is within the linear continuum theory of a uniaxial 
body with a mobile director n that characterises the anisotropy 
directionfl. The free energy density, F, is: 



i§(Tr[!]) 2 - 



±KTr[e](n-e-n) 



Wn-e-n) 



where e = | — jTrJe] 8 is the traceless part of the linear sym- 
metric strain e j; - = j(djUi + djUj), and Q. = ^curl u is the anti- 
symmetric part (the body rotation), u being the displacement 
field. 

The latter terms are those of the relative rotation (ft — uj) x 
n coupling: for small rotations, the director variation corre- 
sponds to a rotation uj = n x Sn w or Sn^ = wxn, see fig. |l[ 
Similarly, a small rotation ft of the matrix causes a vector v in 
the body to suffer the change 8vfj = ft X v. The net rotation 

Sn = uj x n 




FIG. 1: Rotation of n about an axis u>, and body rotations ft rotating 
vectors v. Only Q.2 and fli (not shown) have any effect on vectors 
parallel to n. 

of n with respect to the matrix inflicted by the relative rotation 
ft — uj of the matrix (ft) relative to the changing director (uj) 
accordingly gives a relative change in n of: 

on = (ft — oj) x n . 

The directors n and no before and after application of strain 
are not distinguished in the first terms of equation ([[]) for small 
director rotation, but clearly must be in the relative rotation 
coupling terms in F, 

The tensor e is the strain before it has been made traceless, 



-jjUi [n x e x n] 2 + ^([n x e] -n) 2 
-±Dl[(n-u>) xn] 2 + ±D 2 n-|- [(n-w)xn], 



(1) 



that is, it still has volume changes Tr |^|J in it. Since rubber is 

a soft material, its deformations are at constant volume and we 
neglect the first two terms in equation ([[]) since they involve 



volume change. The geometry of the remaining strains is vital 
and shown in fig. |[ 
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FIG. 2: The elements of strain in a uniaxially anisotropic medium. 
Dashed lines show the undistorted state, whereas the strained ele- 
ments are shown by shaded areas. They divide into stretches E nll 
along n, stretches £,„„, and and distortions e mj (. in the plane per- 
pendicular to n, and distortions e„ m and £„£ encompassing n and the 
perpendicular plane. 

The de Gennes relative rotation couplings are unique to 
nematic networks because they require an independent, rota- 
tional degree of freedom, i.e., the rotations u> of the director 
n. Its motion in the medium is coupled to the body rotation 
ft. Hence, the elastic energy unusually involves antisymmet- 
ric components of shear strain. In section VI , we accordingly 



show how to extract the rotational component from any fi- 
nite shear. The first coupling, D\, purely resists the rotation 
of the director relative to a rotating, undeformed matrix. The 
£>2 term couples n to the symmetric part of shear in the plane 
that involves n (e.g. s ZJ if no = «z in equilibrium). These 
are the shears E nm and e„k in fig. Q Since infinitesimal, in- 
compressible symmetric shear is equivalent to stretch along 
one diagonal and compression along another, it is reasonable 
that prolate molecules will reduce the cost of distortion by 
rotating their ordering direction to being as much as possi- 
ble along the elongation diagonal, depicted in fig. |[ Oblate 
elastomers rotate their director toward the compression diag- 
onal to achieve the appropriate elastic accommodation. The 
spheroid represents the anisotropic shape distribution of the 
crosslinked polymers from which the network is composed. 

A theorem of Golubovic and Lubensky (GL) [|[| shows that 
on symmetry grounds any solid with an internal degree of 
freedom which is also capable of reaching an isotropic ref- 
erence state must be invariant under the double set of rota- 
tions of both the reference and target states when considering 
elastic deformations. This has the remarkable consequence 
that some elastic deformations must be soft. This was discov- 
ered independently and in a superficially different form when 
studying at finite deformations the elastic response of nematic 
elastomers . We sketch how this can be understood for elas- 
tomers in order to explain the complex soft modes we shall 





FIG. 3: Symmetric shear induces director rotation toward the elon- 
gation diagonal. The respective elongations and compressions along 
the diagonals are shown. The chain shape distribution also rotates. 
When mechanical shape changes of a network accommodate the 
rotations of the chain distribution without distortions, such shape 
changes can take place at minimal energy cost. 



later discuss. 

An effective shear modulus arises for an imposed elastic 
strain, when the nematic director is free to evolve optimally. 
The elastic modulus /J2 is reduced; by the symmetry argument 
it must vanish to give no overall energy cost. If the local ro- 
tational component ft of the deformation of the elastic ma- 
trix and the rotation of the director are coaxial, i.e. f2 and 
uj are parallel and uniform, the argument is simple. Min- 
imising the relative-rotation part of the energy density ([!]), 
ADi [(ft — u>) x n] 2 + jL^n-e - [(ft — uj) x n], one obtains the 
optimal relative rotation for a given shear strain e: 

Do £>2 
[(n-w)xn] = -— (n-|)or (fl-v) = -— [n-|xn] . 

(2) 

We now substitute this back into the energy density and obtain 
for the rotation-strain terms: 

iDi[(«-w) xn] 2 + AD 2 n-e- [(fi-w)xn] (3) 



= — n ■ e x nr = — (e xz - 

Wi = 1 8Di v - 



The last expression is written in the specific coordinate frame 
where the initial director no is parallel to z-axis. We can now 
unite this expression with the rest of the elastic energy density, 
equation ([j]) and thus obtain the effective rubber-elastic energy 
which depends only on strains; the director does not appear. 
For instance, in the specific coordinates of fig. 0: 



-2e. l> 2 + e YV 2 ) 



■£V7 



(4) 



The modulus jAi is renormalised to /J2 — D§/(4Z)i) which 
by the GL theorem must be zero, thus establishing a relation 
between the constants /J2, D\, and D^. The molecular model, 
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required below for finite deformations, produces linear con- 
tinuum limiting values which also give: 



D z 



(5) 



Olmsted []6j] first proposed this continuum mechanism be- 
hind the Golubovic-Lubensky theorem: shape depends on the 
orientation of an internal (nematic) degree of freedom, the ro- 
tation of which causes a natural shape change at zero cost 
for suitable solids. A general discussion of the GL argument 
and its extension to semi-softness and thresholds to rotation is 
given in 0]. The experimental evidence for mechanically soft 
distortions is also discussed there. 

One can picture the continuous rotation generating mechan- 
ical distortion at low energy cost, see Fig. |} The natural long 
axis of the body, as defined by the principal axis of molecular 
shape, rotates in an attempt to follow the apparent extension 
axis. To the extent that macroscopic shape change can thus 
be imitated, there is no real accompanying distortion of poly- 
mer shape. We quantify this picture below when considering 
non-linear elasticity theory. 



III. FINITE SOFT ELASTICITY 

A simple extension of classical rubber elasticity theory to 
nematic elastomers gives for the free energy density: 



J F = i i uTr(| -^ T -r 1 -^ 



(6) 



where /j is the rubber elastic shear modulus in the isotropic 
phase and Xjj = dRj/dRoj is the homogeneous deformation 
(gradient). If we take A, to be A, = 8 + u, then the symmetric 
part of u is identical to £ of (jl|) in the infinitesimal limit. 

The chains are no longer characterised by spherical Gaus- 
sian distributions as in the classical case, but in general by 
anisotropic distributions. Thus the effective step length ten- 
sors initially and £ currently after a distortion A are prolate 
(or oblate) spheroids defining the second moments that char- 
acterise the Gaussian distribution of chain spans R in the the 
network, that is (RjRj) = l^ijL where L is the arc length of 
a chain. Thus a measure of the mean size of a chain is £ 1 / 2 
where we shall soon define the roots of tensors more care- 
fully. The tensor £ has one principal value £ii along n and £± 
perpendicular to n; thus 
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where r — £\\/£±- The extracted £j_ factor from £^ cancels 

with the l/£± factor extracted from £~ 1 when both tensors ap- 
pear together in the Trace formula (Q) and we can henceforth 
just consider the £ tensors in their reduced form that simply 
depends on the intrinsic anisotropy r. Anisotropy varies be- 
tween r = 1.1 — > 60. In this model of rubber elasticity the 
spontaneous elongation, A m , on going from the isotropic to 
the nematic phase turns out to be A m = r 1 / 3 and is thus a direct 



measure of r. Indeed spontaneous elongations in the range of 
A m ~ 1 .03 — 4.00 are observed. 

Now distortions A are no longer small, but must continue 
to respect volume conservation, Det(A) = 1 in the non-linear 
regime. The directors no and n, of the initial and current ne- 
matic states, may be greatly rotated from each other. 

The remainder of this paper is concerned with explaining 
the character of the soft modes within finite elasticity theory 
and extending this picture to that of soft modes in biaxial ne- 
matic elastomers. We explore the role of rotations and the 
connections between the isolation of rotational components 
of strain at finite amplitude (the spherical decomposition the- 
orem) and the relation of the roots of tensors to this picture. 



IV. UNIAXIAL SOFTNESS 



Consider the strain I 



\ = £ l l 2 -W 



./-1/2 



(7) 



where W ^ is an arbitrary rotation by an angle a. There are two 
continuous degrees of freedom describing the rotation con- 
necting n and no and three degrees of freedom for the rotation 
W . Hence, the strain A is described by five continuous de- 
grees of freedom and thus represents a large set of deforma- 
tions. 

If we insert such a strain into the Trace formula (Q), as well 
as its transpose A T (equivalent to £q 1//2 • ■ £ 1 / 2 since the £ 
are symmetric) we obtain: 

= ±,uTr(§) = |u. (8) 

Cancelling the middle section, £}l 2 ■ £T X ■ £}l 2 = 8, allows the 
W-terms to meet and gives W T ■ W_ = 8. Likewise disposing of 

the 1L terms, one obtains the final expression F e \ — j/u. This 
is identical to the free energy of an undistorted network. The 
non-trivial set of distortions A of this particular form, equation 
(Q), have not raised the energy of nematic elastomer. From the 
results Det(A-a) = Det(AJDetQ|), Det(A _1 ) = (Det(A))- 1 
and Det(W) = 1 for all rotations W_, one can see that Det(A) = 

Det^ 1 / 2 -^ 'lo" 1/2 ) = 1. *at is these softmodes are volume - 
preserving. 

We saw pictorially in fig. || that, on applying a stretch per- 
pendicular to the initial director, rotation of the chain distribu- 
tion is accommodated by the very elongation we have applied, 
together with a shear. Two remarkable consequences immedi- 
ately follow from nematic elastomer response via rotation: 

Impose an x-extension A xv . All the accompanying distor- 
tions must be in the plane of rotation, that is a transverse con- 
traction A ;; and shears X xz and X zx accommodate the rotation 
of the distribution. No shears perpendicular to this plane, that 
is involving y-direction (A vy ,A vt , .. . etc.) are needed. For a 
classical isotropic elastomer A vv = X zz =1/ is demanded 
by incompressibility, whereas in soft elasticity there is no 
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1/2 



FIG. 4: Chain shape distribution before and after rotations. The ex- 
tent of soft extensions (k xx ) perpendicular to the initial director is set 
by the anisotropy of the molecular shape distribution. The macro- 
scopic dimensions are shown changing affinely with the distribution, 

for instance the x dimension changing from \JT~Z to . 



shrinkage in the y direction (A. vy = 
son ratio is zero. 



1) and the appropriate Pois- 



Fig. |] shows the initial and final states of which the 
shear and partial rotation of fig. [3] is an intermediate step. 
Softness must come to an end when the rotation is com- 
plete and the z dimension has diminished in the proportion 

A,,-. = , and the x dimension extended in the proportion 



'Ka = /£_!_■ The original sizes ^JT^ and \fl~\_ have trans- 
formed to \fl\_ and ^JT^ respectively. Thus softness would 

cease and director rotation be complete at X xx — r 1 / 2 = Am . 
The strain A. m = (£\\/£±) 1 ^ is the spontaneous extension suf- 
fered on cooling to the nematic phase. Likewise one can 
imagine from an oblique form of figure |] that if the initial 
director no (long axis of the shape ellipsoid) is not at 90° to 
the imposed strain, then rotation and softness is complete at a 
smaller X xx < X m . 



We give as a concrete example the set of soft distortions 

1/2 —1/2 

— e i . o i simplified by the absence of the arbitrary 
They are simply characterised (paramet- 
rically) by the angle 8 by which is rotated to |_, that is by 



hoft ~ ie ' *P 
rotation matrix W 



which no is rotated to n. Putting in the dyadic forms for |, 
and I. 



-1/2 



in 



1/2 



into X , gives: 

=soft b 



=sott 



= (8+(V7-l)nn T ).(5+(^ 

= §+(l/Vr-l)nono+(Vr- 
+ (n.n )(2-V7-l/V7)nn r 



■l)nono) 



l)nn n 



(9) 



If no is along z and is rotated by 8 toward x, it becomes n = 
z cos 8 + x sin 8. We can write down a particular representation 



of X . (using the notation s = sin 8 and c = cos 8): 

=soft ° 

X oft = zz T (l-(l--J=> 2 )+xx T (l + (^-l> 2 ) 
T +xz T (l -=)sc + zx T (yfr - l)sc 



-yy 



'l + (V^-l> 2 (\-\/y/r)sc 

1 I . (10) 

(Vr-l)ae \ - {\ - \ / J?)s 2 

The soft modes are neither simple nor pure shear, but a mix- 
ture of the two. Accordingly, the soft modes have an element 
of body rotation in addition to elongations, compressions and 
pure shears. The importance of body rotations is discussed be- 
fore and after (Q). The degree of body rotation is given in sec- 
tion pi Note that the extensional and compressional strains 
(Xxx — 1) and (1 — X zz ) are both proportional to sin 2 8. Thus 
the infinitesimal strains u zz and u xx , at small rotations 8, are 
proportional to 8 2 . By contrast X xz and X zx are proportional 
to sin 8 cos 8 and hence the infinitesimals u xz and u zx are pro- 
portional to 8 - a lower order than u xx and u zz . There is no 
relaxation along y, perpendicular to the plane of rotation n, 
that is Xyy = 1. Also note that in the isotropic limit (r = 1) this 
particular strain X — 8 while the general soft deformation 
^) reduces to the null strain, W , a simple body rotation. Both 
are the trivial cases evidently preserving the elastic energy at 
its minimum. The soft modes become non-trivial deforma- 
tions when the material becomes a nematic elastomer. 

The soft modes start at no strain, X — 8, and as the director 
8 rotates from to Jt/2, they eventually end at 



s/r 
1 
l/y/r 



that is an extension X xx — y/r and a transverse contraction 
Xzz = 1 l\[r. 

The director rotation is taken up by a shape change so that 
there is no entropically expensive deformation of the chain 
distribution as when a conventional elastomer deforms. The 
shape tensor's anisotropy r — t\\/t± characterises the ratio of 
the mean square size along the director to that perpendicular 
to the director. The square root of this ratio, yfr, gives the 
characteristic ratio of average dimensions of chains in the net- 
work. During a soft deformation, the solid must change shape 
such that the rotating ellipsoid i}/ 2 , characterising the physi- 
cal dimensions of the distribution of chains, is accommodated 
without distortion. The ellipsoid is R-^ _1 R = 1, or in a 
principal frame, x 2 /r + z 2 = 1, that is in section the ellipse has 
semi-major axes of %fr and 1. Fig. || illustrates soft deforma- 
tions of a nematic elastomer with chains of anisotropy r = 2. 
The distortions are parameterised by the director rotation, 8, 
which ranges between and 7T./2. The chain shape tensor, 
£ 1//2 , rotates without distortion, just fitting into the solid into 
which it embedded. The soft deformations of fig. |]and equa- 
tions (^|) and ( jlO[ ) of the body are in general non-symmetric. 
Below we decompose these into pure shears plus rotations. 



1 /2 

One can think of X . = £J 

=S0ft =D 



1/2 



as converting the initial 
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FIG. 5: Soft deformations of a nematic elastomer with anisotropy 
;• = 2. The deformations correspond to director rotations of = 0, 
7t/6, Jt/4, 7t/3, 571/12 and 7t/2 which parametrically generate the 
distortions as discussed above. Note that the prolate spheroid char- 
acterising the distribution of chains, embedded in the distorting solid 
(shaded and shown in section), when rotated by 9 can be accommo- 
dated without distortion. The reference, undeformed body is shown 
in outline; the original distribution shown dashed. 

ellipsoid of fig. |]to a sphere by the action of the inverse 1/<2 
and then recreating an ellipsoid at an angle with the action 

of £^ 2 . This is precisely the scheme of [JsJ ^] who consider 
an isotropic reference state (the intermediate created after the 

— 1 12 

action of In ). We discuss this rotational invariance again 
below. 




FIG. 6: The shape ellipsoid for a biaxially nematic polymer. The 
section perpendicular to n is not circular, but has semi-axes of length 
yli p/2. Rotations (]) about n = no generate distortions in the xy 
plane that softly accommodate the non-circular shape as it rotates. 

With fixed n, in the biaxial case one can rotate the 
anisotropic transverse section mk of the square root of the 
shape tensor, £ 1 / 2 , and accommodate it without distortion by 
inducing shape changes in the mk plane. These distortions are 
soft for exactly the same reasons as in the uniaxial case when 
n suffered rotations in the nk or nm planes, see figs. |5|. Given 
that there are (now in the Lab. frame) soft modes as well 
as the X zx and X zy forms (which are now no longer equivalent 
modes), the order of softness has become very much greater. 



V. BIAXIAL SOFTNESS 



A. Explicit examples of soft modes arising from biaxiality 



Biaxial nematic phases are rare. They have been found 
(Finkelmann et al [ |To[ , |llj) in nematic polymers since one 
has more complex molecular structural possibilities. In prin- 
ciple such polymers could be used to make biaxially nematic 
elastomers. They would have rich mechanical properties. 

The shape tensor of a biaxial polymer is 



'l\ o o 
I = ( < 2 
h 



Now the step lengths in the two directions perpendicular to 
n, l\ and 1%, are distinguished. The tensor is reduced by 
taking out a factor of the mean perpendicular step length, 

€x = i(fi+* 2 ),togive: 



I = 



rnn 



l+p/2 
l-p/2 
Or 



(1 



mm 



(l-^)kk T 



= 8- 



l)nn T + ^mm T 



?kk T 



where the axes of the ellipsoid are n, m and k (with the latter 
a third perpendicular axis given by k = n x m), see fig. dq). 



The biaxiality is p = 2 



e !+ e 2 



The situation is exactly parallel to that of the soft shears (|9|) 
and (|Io|). Now the soft modes are given by ^J/ 2 ^ where (]) 
is instead the angle of rotation about the principal director n. 
The tensor £ l J 2 is 



/1/2 . 
SO 



'y/l+p/2 N 

y/l-p/2 




(ID 



which will then be rotated by (j). By analogy with fig. ^, the 
anisotropy in the xy part of the matrix which is to be rotated is 

This measures the anisotropy of ac- 



1+P/2 = 
l-p/2 - 



r± = 



tual (root mean square) chain dimensions in the plane perpen- 
dicular to the director. The resulting soft shears in the xy-plane 
arising from rotations (j) about z = no are: 



=biaxial 



(12) 



More complicated soft shears are possible if these shears are 
combined with those previously found in the uniaxial case. 
For example, in the uniaxial case, section ( |Tv| ) above, subse- 
quent rotation about n after the rotation of n itself had no ef- 
fect since there was symmetry about n. Now that this is lost, 
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such rotations will reorient the perpendicular directors m and 
k and thereby induce additional soft shears. 



B. General biaxial softness 

If, however, the axis of rotation does not coincide with a 
principal axis as in the previous examples, we discover soft 
modes of lower symmetry, which only lead back to an unde- 
formed body after the chain distribution has been rotated by 
271. Consider as an example of reduced symmetry the situ- 
ation when the chain distribution is rotated around a general 
axis in the xy plane at an angle p to the x axis (fig. fy: a rota- 
tion of 7C (fig. |j) transfers the director n (along the z-axis) into 

n, but the unit vectors m and k along the other two principal 
axes end up in a general position in the xy-plane. The same 
final position of the ellipsoid of chain distribution can be ob- 
tained by a rotation around the z-axis by an appropriate angle. 
The corresponding body deformation leaves all z-coordinates 
invariant, but transforms the the xy-coordinates as if they have 
suffered a body rotation of 2p about z. 

If we were to choose the axis of rotation in an arbitrary way, 
we would not observe this feature of invariant z-components 
after a rotation of jc. Only a full rotation of 27C would lead 
back into a state of higher symmetry, in this case, of course, 
the identity. 

By rotating the ellipsoid of chain distribution, we create soft 
modes, which are closely related to the underlying symmetry, 
which is the crystal class of the biaxial ellipsoid, the dipyra- 
midal orthorhombic class. 

Comparing with equation (0), we see that so far we have 
set the rotation W_ to be the identity. If, however, we were to 
allow a general rotation, we would find an even larger set of 
deformations, in fact described by six degrees of freedom: the 
additional rotations and the rotations connecting the prin- 
cipal frames of the initial and final ellipsoids each introduce 
three degrees of freedom. 



VI. ROTATIONS, SYMMETRIC STRAINS AND THE 
ROOTS OF TENSORS 

We have seen in equations (|l])-(0) yielding softness that ro- 
tations are vital to understanding soft elasticity. The cartoon 
of fig. 5 shows that soft shears are not symmetric and hence 
there must be a component of body rotation present. Here 
we explore the role of rotations in symmetry requirements. 
We then explicitly extract the rotational component from gen- 
eral distortions and from soft deformations in particular. The 
requirements for extracting the rotational components of the 
strain tensor are intimately related to finding the square roots 
of tensors which we discuss here since they are used in con- 
structing soft deformations. 

Any general matrix M can be broken down into products 
S ■ O or O ■ 5' of symmetric matrices S and 5' and orthogonal 
matrices O and 0', which can be restricted to proper rotations. 
In effect, one has a pure shear preceded or followed by a body 
rotation (the Polar Decomposition Theorem Jl2[|). 









FIG. 7: Rotation of the chain distribution around an axis (arrow) in 
the xy-plane at an angle of 7t/3 with the x-axis. The figure shows 
intermediate states of rotation at an angle of 0, 7t, i7l,^Jt,^7t,^7l and 
|7t. Fig. (^) shows the final stage of a rotation by an angle of 71. The 
fine outline shows the initial body shape, whereas the highlighted 
outline illustrates the body deformation as the the ellipsoid of the 
internal chain distribution rotates. Here, we have chosen p — 6/5 
and r = 200/45, or, equivalently, the principal axes of the ellipsoid 
to be 6, 3 and 10 respectively. 



Note that in the infinitesimal case, we have 5 = J_+ 55 and 
= L + §Q ? hence 8J = 8J T and §g = . Analogous 

statements hold for the matrices 5' and (9'. This shows that 
the rotation matrices O and & are related to the antisymmetric 
part of the deformation in the infinitesimal limit. 

Applied to the deformation A., we write: 



X = X L -V or U-X* 



(13) 



where the rotations are denoted by V or U_ depending upon 
whether they act on the reference or target spaces respectively 
of the deformation. The form of the accompanying symmetric 
deformations will depend on the order; they are denoted by A. R 
and A. L respectively and yield the Cauchy-Green tensors C and 
B respectively: C = )J • X and B = X ■ )J . 

We summarise classical elasticity theory to highlight the 
differences with nematic rubber elasticity. The deformation 
(gradient) is: 



Xij = dRi/dRoj 



(14) 
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FIG. 8: Rotation of the chain distribution by n around an axis in 
the xy-plane at an angle of Jl/3 with the x-axis. The figure shows 
two views of the same situation: the final ellipsoid together with the 
outlines of the bodies which indicate the corresponding soft mode 
deformation. The outline of the deformed body is highlighted. The 
parameters are the same as in fig. (Q). 

If the target space (St) deforms under rotations represented by 
the matrix U_, as R' = U_ ■ R, and the reference space (Sr) de- 
forms under rotations V_ as R = V_ ■ Ro, then the deformation 
tensor deforms as 

Kj = U ik dR k /dR Q iV7j (15) 

Thus X records the character of both the target and refer- 
ence states, a property that will be essential in non-ideal ne- 
matic elastomers where an isotropic reference state cannot be 
reached. The connection with both spaces is quite different in 
character from that of the Cauchy tensors. Thus the combina- 
tion 

Q = Xj -h = Yj -A' T -g-g T -U'L 

= yj -ti T -X-y = Vj -C^-V^ (17) 

is manifestly invariant under body rotations U_ of the final 
(target) space St and transforms as a second rank tensor in 
Sr. Since isotropic systems are invariant to rotations V of Sr, 
the system's final energy must be invariant to rotations of Sr. 
Thus F is a function of the rotational (in Sr) invariants of C 
and is assured by the above of being invariant under rotations 
of St- As an example, the isotropic rubber elastic free energy, 
setting |° = £ ~ 1 = 8 in equation (||) is 

F = i^Tr (q) = |/iTr ■ g ■ v) 

= ±A<Tr (g ■ Y_ ■ Z T ) = jA'Tr (g) ( 1 8) 

(by cyclical properties of the trace). Likewise, B — X- Xj , is 
invariant to rotations of the reference state and transforms like 
a second rank tensor in the target state. 

Thus we see that "objectivity", frame indifference, is built 
into classical elasticity theory from the outset. Nematic elas- 
tomers are much more subtle. In continuum theory we have 



seen that rotations and symmetric shears enter, both separately 
(through the de Gennes D\ term) and coupled (through the Dj_ 
term). In finite elasticity, we see in the Trace formula that the 
initial and final orientations of the solid and its directors en- 
ters via the tensors ^° and t~ 1 and there are not combinations 
like B_ and C which eliminate rotations. For instance, inserting 

U_ ■ A. R into the trace result, one obtains: 

F = 2MTr(| ^ RT -[£ T -r 1 -£]4 R ) (19) 

where the [. . . ] have been inserted to emphasise that a body 
rotation of St effectively adds to the rotation of n. That is, 
a new £f evolves: t[ ' = Uj£T l U_. (The additional rotation 
is not necessarily coaxial with that which took no to n.) Not 
unexpectedly, we see the effect of U_ compounded with the ro- 
tations implicit in £T X since they both live in the target space. 

Other approaches have been taken [|J U to soft elasticity 
which apparently circumvent the necessity to follow orienta- 
tions in both spaces, and restore objectivity. One can measure 
all deformations A, from an isotropic reference state, that is 
there is encoded into the X first a spontaneous deformation 
to current conditions of temperature, then a deformation im- 
posed with respect to this intermediate state. Under these con- 
ditions the free energy must automatically be invariant under 
operations of V_, since the reference state is isotropic and with- 
out a director no to keep track of. Difficulties arise then when 
nematic elastomers are only semi-soft, that is, they do not de- 
form entirely at constant free energy, because they do not have 
a high temperature isotropic reference state. However they do 
suffer director rotation, very low energy trajectories in X space 
and a lack of X yy relaxation perpendicular to the plane of n's 
rotation JT^]. In these cases, the complete cancellation render- 
ing /jf = m — — > fails, but nevertheless deformations are 
qualitatively soft and one has to keep track of both directions 
no and n as strain evolves. 



A. Square roots of tensors and the Polar Decomposition 
Theorem 

We quote the classical conditions jl^] for the square roots of 
tensors, that arise too in the all-important condition for polar 
decomposition that is required for nematic elastomers. 

If A is non-singular, with fj distinct eigenvalues and with v 
Jordan blocks, then A has > 2 M and < 2 V non-singular square 
roots. At least one of these roots is a polynomial in A. 

The proof of the Polar Decomposition Theorem (PDT) also 
offers a practical algorithm for decomposition: take a non- 
singular X and construct the manifestly symmetric XjX. Take 
the square root, G, of A. T A. that is a polynomial in A. T A.. Since 

}7X is symmetric, then so is G = ^JxJX. Then define Q = 
G X. Clearly Q T Q — 8, that is Q represents rotations. From 
this one recovers X= GQ. 
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B. Symmetric strain and body rotation of soft modes 

We investigate the special case of uniaxial soft modes in the 
xz plane. They are achieved by rotating the chain distribution 
around the the y axis. The resulting deformation keeps the y- 
components of the body constant. Hence these soft modes are 
effectively 2x2 and can be easily related to equation ( |ic| ) and 
fig- 1 

An explicit example of the PDT for decomposing a general 
xz-distortion A into a combination of symmetric distortion, A, , 
followed by a body rotation, U_, about the y-axis by an angle 
a is: ~~ 



=soft 



Xxx 8 
5' JU 



c s 

— S C 



(with s = sin a and c = cos a). One can confirm that the rota- 
tion is: 



tana = (8- 8')/ (A* 
and thus for instance for sin a: 



(20) 



sin a = — (8 - 8') with A 2 
A 



(A xr + AJ z + (8-8') 



l\2 



Of course there is no body rotation for A symmetric, that is 
8 = 8'. One can also confirm that the symmetric shear tensor 
is : 



A R = 



1A„ 
A V 



fe+M-8'(8-8') Xxxb+S'Xx 
A. xr 5+8'?i = \ zz (k zz +\ xx )+5(8-8') 



(21) 



The results ( ^0| ) and ( |2l| ) break the soft modes down into a 
symmetric shear A R followed by a body rotation U_ through an 

angle a about the y-axis. Thus the soft mode is A soft = U_ ■ A R . 
We continue to parameterise them with the director rotation 8. 
The body rotation is through an angle a given by 



tana 



sin6cos6(V?- I) 2 
2V^+sin 2 9(V?-l) 2 

tan9(/F- l) 2 
2^/7 + tan 2 8(> + l) ' 



(22) 



For small shear (small director rotation 0), the component 
of body rotation is small too and proportional to 9: 



oc«8 



(V^-i) 2 



(23) 



For large rotations, as Jt/2, the rotation a vanishes as we 
have seen in fig. |J Thus it is only at first that body rotation 
plays a part in accommodating the rotating chains. As the ro- 
tation of chains approaches 7c/2, the body is simply extended 
or compressed along the principal axes of the original chain 
distribution. 

The corresponding symmetric shear strain d, the off- 
diagonal component of A s , is 



1 sin0 cos8(r 
A 



1) 



which, for small distortions is also proportional to 8, bearing 
in mind that A is constant at first order in 8: 



1 80-1) 
A V? ' 



(24) 



As we have seen in fig. |[ it also vanishes at 8 = 7t/2, where 
no further accommodation of the shape tensor by rotation is 
possible. In fig. ^, we show the rotation and off-diagonal 
element of the symmetric shear occurring during soft defor- 
mations between 8 = and K. Notice that both quantities are 




FIG. 9: Off-diagonal element d of the symmetric shear A R , and an- 
gle a of the rotation U_ (dotted line) plotted against director rotation 
6 for the soft deformations of fig. n for a nematic elastomer with 
anisotropy r = 20. 

linear in 8 for small distortions. As the director of the chain 
distributions is rotated beyond k/2, the solid body rotation 
a and the off-diagonal element of the shear d become nega- 
tive. At a director rotation of Jt, the original body shape is 
recovered, as it should be in an uniaxial nematic system under 
reflection of the director n — ► n. 

Similarly, the behaviour of the pure shear is instructive: for 
this purpose, we plot in fig. |l0| the two eigenvalues of A , 
which serve as a good measure of the net presence of shear. 
Note that, due to incompressibility, the product of the two 




FIG. 10: The two eigenvalues of the pure shear A for a nematic 
elastomer with anisotropy r = 20. For a director rotation of 7t/2, 
the soft mode deformation becomes a pure shear with extension and 
compression ratio yfr or \j \fr respectively. 

eigenvalues is constant to 1 . Both eigenvalues are linear for 
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small 8, saturate at values y/r and 1 / ^fr respectively for a di- 
rector rotation of 7t/2, where the soft mode becomes a simple 
compression and extension along the principal axes. At a ro- 
tation of Jt, we recover the eigenvalues of the identity matrix. 

^From equation (|l0|), we see that the compression and ex- 
tension in the x and y directions are quadratic in 8 for small 
director rotations. In other words, u zz and u xx too are propor- 
tional to 8 2 oe a 2 d 2 (by equations @ and (§|)). 

VII. CONCLUSIONS 

Soft deformations of nematic elastomers result from the ro- 
tation of the anisotropic chain shape distribution without dis- 
tortion and therefore without rubber elastic free energy cost. 
This is by contrast to classical rubber, where distorsion of the 
distribution lowers entropy and raises the free energy. Ex- 
plicit forms of these soft modes are derived in the case of 
biaxial nematic elastomers using the Olmsted method of the 
square roots of tensors developed for the uniaxial case. For 



both cases we show geometrically what these deformations 
look like. They correspond to rotations of prolate distribu- 
tions being accommodated by elastic strains of the body they 
are inscribed into. 

The elasticity of nematic elastomers depends, unlike in 
classical elastomers, on body rotations (since these can be 
with respect to the underlying nematic director). At finite 
strains it is important to isolate the rotational component of 
the generally non-symmetric deformations. We make contact 
with the Polar Decomposition Theorem in this context. In do- 
ing so we discuss the roots of tensors that are employed in 
finding the manifold of soft deformation tensors. 
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